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ABSTRACT 


In this thesis some special cases of the rapacity 
expansion in the flow network are considered. Specifically 
the cases considered ares (i) From the given set of arcs 
which can be added to an existing network, select an arc such 
that it maximizes the increase in flow from source node to sink 
node, (ii) From the same set of arcs select arcs with minimum 
total capacity such that their addition to an existing network 
increases the flow from source node to sink node by a desired 
value . 

For the first problem an implicit enumeration algorithm 
fdr which identifies the optimal arc is developed. For the 
second problem two heuristics and one exact method have been 
suggested. The exact method is an implicit enumeration scheme. 

The heuristics utilize the lower and upper bounds on capacity 
utilization of the arcs to be added. 

The computational performance of the enumerative algorithm 
(for single are addition problem), implicit enumeration scheme 
and one heuristic method is tested on randomly generated networks. 



CHAPTER I 


INTRO DUCTION, 


1.1 INTRODUCTIONS 

Physical distribution systems such as electricity or gas 
supply, road or rail networks, etc. can be easily modelled as 
networks of flow. Generally these systems are designed optimally 
with respect to the existing requirement. However with changes 
over time additions or contractions of these flow networks may 
be required. . 

After some time it may be desired to increase the flow 
in this network to a higher level. This can be done by either 
increasing the existing arc capacities or by introducing new 
arcs (or both). However, increase of this capacity will be at 
additional cost and as these costs are substantial, arcs with 
minimum cost would be the most desirable ones to be introduced 
in the network. This problem is very well known problem and 
considerable amount of literature is available on the capacity 
expansion of network problem. 

However, there are cases of the problems where the capacity 
expansion may be done under restricted conditions. One simple 
case is when the costs of the arcs added are linear function of 
the capacity of the arcs. In such cases the problem of adding 
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arcs with minimum cost will be equivalent to the problem of 
adding arcs with total minimum capacity. 

Similarly, we may be interested in finding the arc from 
the list of specified arcs whose addition will increase the flow 
to maximum. Such an arc can be viewed as most valuable arc. 

In this thesis, we will concentrate on the above men- 
tioned two problems. In the next few sections, we will review 
the literature in the related fields. 

1.2 LITERATURE REVIEW? 

There is considerable amount of literature available for 
the capacity expansion problem and related areas. For convenience 

we will review literature in the following areas. 

(i) Maximum flow problem. 

(ii) Minimum cardinality min-cut problem. 

(iii) Network sensitivity analysis'. - 

(iv) Capacity expansion problem. 

1.21 Maximum Flow Problem ? 

Maximum flow problem can be briefly described as belowi 

Given a network G(V, A) with vertex set V and arc set A and 

capacity c. . on arc (i,j) find the maximum flow from a specified 
id 

source node s to a specified sink node t such that flow is 
conserved at all the nodes except at source and sink and no arc 
has flow larger than its capacity. It is a well known problem 
and is a special case of linear program. Ford and Fulkerson s [ 3 ] 
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labelling algorithm was the first algorithm proposed for this 
problem, which essentially finds maximum- flow through a series 
of augmenting paths. Algorithm is finite and gives a integer 
flow in the case when arc capacities are integer. 

Dinic 12] has proposed an improvement over this algorithm. 
Dinic' s algorithm involves constructing layered networks with 
respect to existing flow. These layered networks are the union 
of all shortest paths (in terms of number of links iri path) 
between s and t, such that each path has flow augmenting capacity. 

A maximal flow in this network is obtained. The process is repeated 
till the maximum flow is obtained. The complexity of algorithm 
is 0(|v|^), where fV j denotes the cardinality of the vertex set V. 

Karzanov [ 8 ] has further improved the method of finding 
maximal flow in the layered network, using concept of preflows. 

This algorithm has a overall complexity of 0(jv|f). A similar 
but much simpler algorithm with complexity 0(|v|^) is proposed 
by Kumar, Malhotra, Maheshwari [13]. A recent study [6] has 
compared the various algorithms for its computational efficiency 
and it concludes that Karzanov 7 s and Dinic r s algorithm perform 
better than other algorithms on average on randomly generated 
networks . 

1.22 Network Sensitivity Analysis s 

M«*n -t «nr'i»nMV«< -*.>r .nu nc ■ ■ mnn, r mm * 

Problems of arc additions and arc deletions can be viewed 
as problems of sensitivity analysis on the network. 
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A most vital link in a flow network is the arc that 
reduces the flow by maximum amount when it is deleted from 
the network. Wblmer [l6] has developed an algorithm to obtain 
the most vital link. He has also given the following necessary 
condition for a link to be a most vital links 

A necessary condition for an arc (a, b) to be most vital 
link is that for any maximum flow pattern in the network G(V, A) 
the flow in the arc (a, b) is at least as great as the flow 
over every arc in a minimum cut. An improvement on the algorithm 
developed by Wolmer has been made by Lubore, Ratliff, Siclia [12] • 
They identify the set of arcs satisfying the above necessary 
condition as candidate set. Then they have developed a labelling 
scheme similar to that of Ford and Fulkerson s / labelling algorithm 

for obtaining maximum flow between given source and sink. This 

0 

labelling’ scheme systematically searches for a flow path from 
node a to node b which does not include arc (a, b) such that 
flow may be diverted from arc (a, b) over this path where (a, b) 
is the arc from candidate set. If node b is labelled then a 
breakthrough occurs. The occurrence or nonoccurrance of breakthio ugl 

is used to ultimately determine the most vital link. 

% 

The n arcs in a network whose simultaneous removal from a 
connected single commodity flow network results in the greatest 
decrease in the throughput capacity of the remaining system 
between a specified pair of nodes (source and sink) are the n most 
vital arcs of the flow network. 
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Ratliff, Lubore, Scilia [15] have developed a efficient 
method to determine n most vital links of the network. 

More precisely the problem of finding the n most vital 
link can be stated as that of finding a set of n arcs D* in 
G(N, A) such that min-cut separating s and t in G(N, A-D ) is 
less than or equal to min cut separating s and t in G(N, A-D^) 
for all sets D n c: A with { D } = n. 

A modified network procedure is developed first. Then 
the general algorithm which uses this modified network procedure 
is developed. If the modified network procedure fails to obtain 
the solution in one pass then General algorithm is used which 
uses modified network procedure several times to obtain the 
solution. It may be noted the problem of adding n arcs to 
increase max flow can be related with the removal of n arcs from 
the complete graph. 

1 » 23 Mi n Cardinality M in Cu t Probl em; 

The problem to be considered is finding minimum cuts in 
the given network with the additional property that the number of 
arcs in the cut (X, X) or (3?, X) or (X, X) U (X, X) is minimized 
where (X, X) U (X, X) is minimal cut of the given network. 

Mamacher [ 5 ] has given a method for obtaining the min 
cardinality min cut. The theory of lexicographical network flows 
has been employed to develop the algorithm. 
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1 . 24 Network Capacity Expansion Problem s 

Fulkerson [ 4 ] has considered the problem of allocating a 
budget of resources among the arcs of a network for the purpose 
of increasing its flow capacity relative to given sources and 
sinks. The assumption made is that the cost of increasing each 
arc capacity is linear. A labelling algorithm is described that 
permits rapid calculation of optimal allocation for all budgets. 

The above assumption permits direct formulation of the problem 
described above as a linear program. The algorithm produces 
solution to the problem, not only for a fixed budget, but for 
all budgets i.e., the problem is solved parametrically. The 
algorithm uses a variant of the labelling procedure previously 
developed to solve the maximum flow problems. 

The long range planning model of transportation system to 
assist planners in assessing the impact of various levels of 
service in transportation network has been developed by Le-blanc[llJ 
In particular rail road network distribution are considered, the 
emphasis being placed on the problem of determining optimal levels 
of service on each arc in the existing network. 

The model is formulated so that improvement to the shipping 
arcs are included as decision variables. The concave transhipment 
problem is extended to the case where co-efficients of the shipp- 
ing cost functions are treated as decision variables. 
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Thus the model determines optimum movement of freight 
between nodes and optimum improvements/degradations to the net- 
work to minimize shipping costs plus maintenance costs on the 
shipping arcs. 

Christof ides and Brook er [ 1 ] have considered the problem 
of increasing the capacity of existing network where new arcs 
of finite capacity can be added at certain cost. The problem 
considered is -Given an existing network, a list of arcs which 
could be added to the network, the arc costs and capacities and 
an available budget, the set of arcs are to be selected such 
that their addition will maximize the maximum flow from source 
s to sink t subject to the budgetary constraints. 

The method describes an efficient tree search algorithm 
using bounds calculated by a dynamic programming procedure which 
are very effective in limiting the solution space explicitly 
searched. Computational results for a number of medium sized 
problems are described and computing times are reasonable. 

Howard and Nemhauser [ 7 ] have considered the following 
problems Given a sequence of predicted demands for N time periods, 
determine the optimal investment decision in each period to mini- 
mize a linear investment cost and a strictly convex cost capacity. 
The relationship between capacity and the investment decision is 
assumed to be linear, but time varying constraints on both the 
individual decisions and on the sum of decisions are considered. 
They have also derived an algorithm for solving this, problem. 



CHAPTER II 


SINGLE ARC ADDITION IN A FLOW NETWORK 

2 . 1 INTRODUCT IONs 

In this chapter we will consider the problem of adding 
a single arc in the flow network G(V, A) with a specified source 
vertex s and specified sink vertex t. We will like to select 
the arc to be added such that the increase in flow between s and 
t is maximized. Such a problem may be meaningful in case where 
the cost of addition is identical for all the arcs and the budget 
permits only addition of a single Karo. 

It is obvious that if such an arc is added between source 
and sink directly, it will result in the maximum increase in the 
flow. However, the choice of arcs to be added may be limited to 
set E and such an arc may not be there. Thus it will be worth- 
while to explore the addition of arc in the case where no such 
s, t arc is permitted to be added, or the capacity of such an arc 
is small. In such case identification of such a link can be 
considered as most valuable link as its addition will increase 
the flow to maximum. 

2.2 DEVELOPMENT OF THE PROBLEMS 

We will consider a network: G(V, A) where V is the node 
set and A is the arc set. Vertex s will denote the source and 
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vertex t the sink vertex. All arcs are undirected. Let c^j 
denote the capacity permitted on the arc (i,j). A cut denoted 
by Y is a partition of the node set into two disjoint parti- 
tions Y and Y such that s e Y and t e Y. Further capacity of 
cut Y is denoted as 

C (Y, V) = E c 

' ieY 10 

jeY 

The flow in this network is governed by the following theorem 2 
given by Ford and Fulkerson [ 3 ] . 

Theorem 2.1 (Maxflow min cut theorem) 3 

Max flow f Min C(Y, V) 

y 

An obvious corollary to above theorem is as follows. 

An arc (i, j) will add to the max flow only ifEit spans 
all the minimum cuts in the network, i.e. it adds to capacity 
of all minimum cuts. 

It may be noted that if(s, t) arc is permitted to be 
added, then it will span all the cuts and hence will be the arc 
to increase the flow to the maximum extent. 

We shall concentrate on the restricted version of the 
problem, where arc to be added is to be selected from the set E 
of arc- and the capacity of each arc is specified. 
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A simple enumerative scheme for solving this problem can 
be constructed by solving max flow problems after adding each 
arc individually and selecting the arc which increases the flow 
to a maximum value. However, the procedure will require 
0( |Ej |vp) effort. Enumeration can be substantially reduced 
by using corollary 2.2. As the selected arc has to span all 
the minimum cuts one can reduce the computation by obtaining 
a minimum cut Y, Y and selecting from the arc set E„ set of 
admissible arcs as E' such that (i, 3) e E'if i e Y, 3 e Y. 

Let f. • be the flow on the arc (i, 3) corresponding to max 

J. J 

flow f. Let denote the node potential of node i i.e. amount 
of flow which can be sent from node s to node i e Y or from node 
t to i i e Y , in the referant network. The referant network is 
constructed by adding two directed arcs corresponding to each 
undirected arc (i, 3) with flow f. The capacity of edge (1,3) 
will be c. . - f . . and of ware (3,1) will be c . . + f. . in the 

A-O 

referant network. Then the maximum flow which can be increased 

by adding an arc (m, n), denoted by p m ^ n = Min (P m , P n , c^}. Also 

the flow will increase to maximum by adding the arc (p, q) which 

is max min (F , P n , c^) . In case E ' contains arcs of the type 
(m»n) e E 

(ij t) for i e Y and (s, 3) for 3 e Y and the arcs are uncapacitated 
the maximum increase in flow will be on one (i, t) or (s, i) as 
the case may be. 
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Further reduction in computation can be obtained by identi- 
fying a minimum cardinality s cut and a minimum cardinality t cut. 

A cut Y is minimum cardinality s cut, if j Y | is minimum among all 
minimum cuts. Similarly a cut Y is minimum cardinality t rut 
if |Y| is minimum among all such cuts. Let such cuts be denoted 
by (S, S) and (T, T) respectively. 

Corolla ry 2. 

An arc will increase the flow only if it is of type (i, j), 
i e S, j e T. Thus E can be further reduced by selecting only 
those arcs which are of the above type. 

Further in the following lemma, wee show that min cardina- 
lity cuts (S, S) and (T, T) are unique. 

Lemma 2.4 s 

Minimum cardinality cut (S, S) (and therefore (T, Tj ) is 

unique . 

Proofs 

«•*£>*<****» «ttUK 

Proof follows from the following theorem, as stated in 
Ford and Fulkerson [ 3 ] . 

Th eorem 2 .5s 

If (X, X) and (Y, Y) are two minimum cuts then (X (1 Y), X'TTy) 
Is also a minimum cut. 

If (S, S), the minimum cardinality cut is not unique, then 
;here exists another minimum £ cordinality^cut (S’, S'). Further 
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|S n S' I < I s j and (S n S', S n s') is also a minimum cut, 
contradicting the fact that S is a minimum cardinality cut. 

Next we describe a procedure to identify a minimum 
cardinality cut. 

2.3 ALGORITHM FOR MINIMUM CARDINALITY CUTS 
Proc edure MINCARN g 

Step is Solve a s, t max flow problem on the given network; 
obtain a min cut (Y, Y), 

St ep 2} Increase the capacities of all the arcs to nc p-j where 
n= }V| (number of nodes in the network). 

S tep 3 » For all nodes i s Y add a arc (i, t) with c^ = 1 

where c^ is the capacity of arc (i, t) or if there 
exists an arc (i, t), increase the capacity of arc (i, t) 
by 1. 

Step 4; Solve a s, t max flow problem. The new min cut (S, S), 
thus obtained will be min cut in the original network 
with the additional property that |o| is minimum. 

Proof, s The capacity of the minimum cut in the modified network 
will be due to two types of arcs. 

(i) The arcs already present in the original network. 

(ii) The new arcs added (of (i, t) type) to modify the 


network. 
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Let the modified network be G ' . Let the min cut in the 
original network G be K and in the modified network G f be K' . 

Let us assume that the min cut in the modified network is not 

a min cut in the original network. 

As | K ’ j >,1, there will be at least one arc from node in 
K' to node t of capacity 1. Also since all arc capacities are 
integers, capacity of the cut K* in G', 

C (K\ G 1 ) > nc (K\ G ) + 1 (2.1) 

where C (K', G ) is the capacity of cut in K ' in G. 

Also C (K, G') < nc (K, G) + n - 1 as |K| < n-1 (2.2) 

where C (K, G') is the capacity of cut K in graph G'. Now as K 1 
is not a minimum cut in G , 

C (K\ G ) > C (K, G ) + 1 (2.3) 

Hence from (2.1) and (2.3), C (K* , G' ) >nC(K f , G ) + 1 

> nc (K, G) + n -f 1 (2.4) 

Thus we have from (2.2) and (2.4), 

C (K* , G r ) > C (K , G') 

which contradicts the fact K 1 is minimum cut in G ' . Therefore K* 
will be minimum cut in G ' only if it is a minimum cut in G . 

Further, C <K\ G’) = nC (K, G) + j K 1 | and as K* 
is minimum cut in G' it will be the cut with minimum | K' | i.e. 
minimum cardinality cut in G. 
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2.4 UPPER AMD LOWER BOUNDS ON ADDITIONAL FLOW BY ADDITION 
OF AN ARC; 

As outlined in Section 2.2 the enumeration scheme for 
identifying the arc to be added to increase the maximum flow, 
will involve solving one max flow problem for each arc in E . To 
reduce the complexity further we will explore the idea of asso- 
ciating a priority order for the addition of these arcs. For 
this we will identify procedures for computing lower and upper 
bound on the additional flow which can be obtained when a arc 
is added. We will identify lower and upper bound on the flow 
permitted by the nodes, for all nodes in minimum cardinality s 
cut, or minimum cardinality t cut and then compute the lower bounds 
on additional flow which can be obtained when specific arc is 
added. 

2 • 41 Identi fication of Lower Bound s o n Node^Pc ^entiali 

Given a maxflow pattern in the existing network correspond- 
ing to minimum cardinality cut (S ? S) the flow potential for each 
node in S can be obtained by solving maxflow problems between s 
and the node i in the network with residual capacities. However, 
this procedure will require solving a maxflow problem for each 
such node and will be time consuming. lower bounds on each node 
potential can be obtained by identifying a good feasible flow 
between s and the node in the network with residual capacities. 

We will describe two procedures which obtain such flows. 
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2.42 Pro cedur LAYERS 

In this procedure we shall first obtain a minimum cardi- 
nality s cut (S, S) and a minimum cardinality t cut (T, T) of 
the network G (V, A) . Then a sequence of non-crossing cuts 
S = S k , S k ^ , ... s’"' = s are obtained in a reduced network where 
all the nodes of S are condensed into a single sink node. 

Further, let V^, V k _^, v c denote the partitions 

of S such that j - U V . It may be noted V Q = {s} and 

V x nV 0 -*• 

Let C(S^) denote the capacity of cut (S , S ) in the 
original network. Then these non-crossing cuts are obtained 
such as to satisfy following conditions s 

(1) C(S°) > c(s' ) > ... > C (S k ) 

(2) The node partitions are such that for any arc (i, j) e a 
if I f s, is in the node partition V^ s then j is in either node 
partition V or in V (or V ). It implies that except from 
the source node all arcs will be cither within same node parti- 
tion or in adjacent node partition. 

To obtain these cuts we iteratively follow following 
condensation procedure. If at the current step the current cut 
is (S 1 , S*) in the network G (V, A) then condense the nodes of 
S* in a single node t i . Remove all the arcs from A which are 
in the cut (S*, S i ). For each node p e S 1 , if there is a arc in 



16 


* * * 

cut (S 1 , S 1 ) add an arc (p„ t 1 } with capacitor c . =■ M (a 

p * 1 

large number >_ 2 c . . ) . 

(i,j)eA 

Solve a max flow problem between s and t 1 . The minimum 
out obtained corresponding to this maxflow will be S^ \ Repeat 
the procedure till S° = s. 


In the above procedure following may bo noted. 

(i) S 1 " 1 c: S 1 

(ii) If there is a node p e which has an arc to a node 

q e S^, j > i then node p 0 S 1 " 1 except when i = s. It follows 
from the fact that this node p will have an arc (p, t 1 ) with 
capacity M in the reduced network corresponding to (S , S> ) , and 
can never be capacitated, this implies that p e S , and hence 
p jZ S 1 " 1 . However, if p = s then s having unlimited capacity 
will belong to S 1_1 

(iii) Partition V ± = S 1 - S 1 - 1 - 


(iv) v i+1 n v i = 0 

(v) From (ii) it follows that for any 

then q e V ± or or V ±+1 except when p 

(vi) Further C (S i_1 ) >. C (S^) as ( S i 
condensed network in which minimum cut S 


arc (p, q) if p £ 

= s. 

t 3 i_1 ) is a cut in the 
was obtained. 


Similarly corresponding to minimum cardinality t cut (T, T) 


shrill obtain a sequence of noncrossing cuts T°, T^, 


we 


T 
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such that T° = T,..., T m = t. These cuts will also have the 
properties similar' to properties of stated above. These cuts 
can be obtained by the procedure outlined above except that t 

will replace s in the procedure. 

. 4 

The resulting network will have the following structure 



2.43 AlgoritlM.,LAYgRs 

Inputs Network G (V ? A) with capacities c^^ of arc (i, j) and 
with specified source node s and sink node t. 

Ste p Os Find the min cardinality s-min cut in G(V, A) using the 

algorithm MINCARN of Section 2.3. Let S be the partition 
set of the cut such that s e S. 
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k-*= — • 0 

X — s , x -fr - - V - s 

Comments X is the current cut partition containing node s. 
Step, is For all (i, j) s.t. i e X, j e X, 

C ij*'~~ ° 

For all i e X s.t. there is an arc (i, o) s.t. j e X, 
add an arc (i, t), 

c^ «*s — M (Comments M is large number). 


Stop 2 s Solve maxflow between s, t and let the minimum cut 
obtained be (S, S) then, 

"v. < — X - S 
k 

X S, X -6 V - S 

IF |XJ = 1 GO TO Step 3, else k<. • k+1, GO TO Step 1. 

Step,J_s For i * 0, k, V.* 

Stepjis Output % V o , V 1# ..., V k . 

Similarly the node partition of the minimum cardinality 
t cut can be obtained by interchanging s and t. STOP 

2 * ^ Lower Bou n ds r on N od _e_ P otential Using r _Pro endure _ LAYER s 


Once the layered network and the node partitions 

V. , v i if . .., v t v °» **•> arc obtained, the lower bounds 

K. w K— X w O 

on the potential of nodes i e S find i e T, can be computed in 
the following fashions 
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For simplicity we shall renumber the nodes such that 

for all arcs if i e V and 3 e V , then i < 3, and arc is 

P P + x 

regardo-d as (i, 3). 

Compute for the each undirected arc (i, 3), the residual 

r 

capacity with respect to maximum flow as follows % 



c . . - f. . if the flow on arc (i, 3) is in 
direction (i to 3). 

c . . + f . . if the flow on the arc (i, 3) is 
-*-o 0 

in direction (3 to i). 


Compute the lower bound JL on the node 3 e V as follows* 

J P 


Let a. 

u 


= I r. . + r . 
. 1C S3 

(it 3) l ieV p -i 


then, 

L = a. + Min {a., r, -} 

1 J i £ v x ij 

a. denotes the amount of flow which can be pushed to node 

3 directly from the nodes in the previous layers. As the layer 

V is constructed by shrinking the nodes in ST , there exists 
P P 

feasible flow such that all edges between V +1 , will be 
saturated, and hence node 3 can receive at least additional 
flow a.. Further from each node i in the same layer it can 

tj 

get flow which is minimum of and r.y Thus are feasible 
additional flows from node s to node 3. 



20 


2.45 Procedure REFERANT Networks 


A simple procedure for constructing a layered network 
can be developed based on Dinics [2 ] method for max fl~w prob- 
lem. In this procedure, first a maximum s, t flow is obtained. 
Next a directed network having two arcs corresponding to each 
arc in the original network is constructed. For convenience - " 
we shall label all undirected arcs as labelled in the direction 
of the flow. If f^j is the flow on arc (i, j) from node i to 
node j and the capacity of this arc is c. then the new network 

-L J 


will have two arcs (i, j) and ( 3 , i) withe capacity (c 


id 


V 


and (c^ •+ f^j) respectively. If for any arc in this network 
the residual capacity is zero, the arc is removed from the network. 


Construct the layers V , V , . ... V of the nodes as 

J o jl i n 

follows s 

V o - < s > 

V p “ hi (1,3). ie Vl* 3 fZ lio Vk } 

Let V , V 1 , ..., be the partitions of the min cardina- 
lity s-min cut, constructed by this procedure. A simple lower 
bound can be obtained by computing the ore disjoint path from 
s to node 3 and then adding the residual capacity of each arc 
disjoint path. The residual capacity of:, arc disjoint path is 
the residual capacity of minimum residual capacity arc in the 
path. 
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Thus £. = £ r where P is the set of arc disjoint 

J peP P 

path between s and j and r denotes the minimum residual capacity 

ir 

(augmentation capacity) of path P. 

We next describe an algorithm to obtain these lower 

bounds. 


2 . 46 Algorithm REFERANT Networks 

Inputs G( V, A) with c . . as arc capacity ofarc (i, j) and 

-4 J 

specified source and sink. Compute maxflow f and 

flow f . . on the arcs using Algorithm MINCARN. Let 
3-d 

(S, S) be the min cardinality s cut. 


St m ep_JLj 


For every arc (i, j ) e A 

If f . - >0 then c— c . i 

3-d ’ 3*d 




+ 



v S . V'- ■— 1 

o J 


A 


*■ 0 


StepL.2» V p ^ — 0 , k * 0 

Repeat for all (i, j) such that i £ V ^ 
If ®i j > 0 and 5 ^ v p-i 

V" V P u 

k <-■ — k+i 
JC — A U { i, j } 
else delete arc (i, j). 

If k =0, m<— -p, GO 10 Step 4. 
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St e&Jtf p«*— p+1, GOTO Step 2. 


Step 4; Output V Q , V ffi , S, A. STOP. 

Similarly algorithm is repeated to obtain partitions 
V^, V° of minimum cardinality t cut T. 

2.47 Procedure RLBONDs 


Comment: This procedure computes the lower bounds on the node 
potential of nodes in S (or in T) . The lower bound on the node 
j is the sum of the residual capacity (flow augmentation capacity) 
of arc disjoint path from s to p 


INPUT: Network G(V, A), Capacity of arc (i, j). 

Step As Using procedures KEFERANT Network obtain V^, V^, 


Step 0 i 


V s > c id-’ *id 

p * — 1 




B .<- ■■ 0, for all j £ V 
1 


id 


0, for all i, j e A 


§£eEJL* Find a node C e v p such that B j = °* 
z « — m (Comments large number) 

GO TO Step 2. 

IF no node j e such that B.. = 0 is available, 
GO TO Stop 7. 

step,z £ q*~~ 3, 1 
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S£§RJ3s Find an arc (k, q) such that k e V a and F =0 

p-i, pq 

2 min {z, r kq } 

F. <e— 1 
kq 

GO TO Step 4. 

IF no such arc is available, GO TO Step 6. 

SLep. kt if~~ 2+1, q — -k 

IF (p < 0), GO TO Step 5, else GO TO Step 3. 

Stex> ,5.8 GO TO step 2. 

Step.. 6.8 F kq ^ — 0, for all (k, q) £ A 

B <r 1 

J 

GO TO Step 1. 

Stepps p«— p+1 

IF p < m, GO TO Step 1, else GO TO Step 8. 

Comments m is the number of layers in the network G ( V, A ) . 

Sj^ep_.gs _ Lower bounds on nodes have been computed output a' .STOP. 


2,48 Upperbound on the Node Potentials 

$ 

Simple upperbound for each node i e S can be computed by 
noticing that the maximum amount of flow a node can receive from 
the source is equal to sum of the capacity of arcs incident on 
this node which have one end in min cut partition S less the 
sum of the capacity of arcs which have one end in min cut 
partition S and are incident on this node i.e. 
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u. = 2 c. . - Z c. . 

1 , . 10 lj 

(i» d) £ B^ns (i,d)sB i ns 

where is the set of arcs which are incident on node i. 

2.49 Upper and Lower Bounds on the Arc si 

mmMn *"*'-*■< "tuirysowauT in ~~r **«i.i**i-*s:-r»* tms utstm , - t . — 1 . r t^M ik — - - -T-, - 

Let E 1 constitute the feasible set of arcs from E whibh 
can be added to network G (V, A) . 

Arc potential of an arc of E 1 is the amount of flow it 
can increase in the network. 


We obtain the lower and upper bound on arc potential for 
every arc of E' in the following fashion. 

Lower bound on arc potential of arc (i, j) e E* is given 
by the following equation. 


min { » 


P. 0 . 


c. .} 
1J J 


where P . is the lower bound on node potential of node j and 
''J 

c ij is the ca P acit y of arc (i* d)* 

Upper bound on arc potential of arc (i, j) e E 1 is given 
by the following equations 


v io “ Min { u i’ u j’ c ia> 

v^iere u. is the upper bound on node potential of node J. 
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2.50 Proc edure for Determining the Arc Causin g Maximum 
Increase in Flow: 


In this section we give a method for selecting an arc from 
the given set of arcs E , such that when it is added to the 
network G(V, a) the increase in flow from s to t is maximum. 

In this enumerative scheme arcs in E , will be arranged 
by a given priority rule in a set Q and will be selected from 
Q one by one. The four priority orders we suggest are as 
follows: 


(a) Q=E / i.e. arcs in Q are arranged in the same order as 


that in E 1 . 

(b) a is the set of all arcs e E' which are ordered in the 
non increasing order of the lower bounds on arc potential. 

(c) Q is the set of all arcs e E / which are ordered in the non 
increasing order of the upper bounds on arc potential. 

(d) Q is the set of all arcs in E' which are ordered in the 
non increasing order of the capacities of arcs. 


At a general step an arc from Q is selected and its upper 
bound an arc potential is tested with the incumbent solution. 
Depending upon whether the upper bound is greater than the incum- 
bent solution or not the arc is added to the network or pruned 
from further consideration. The incumbent solution is dul> 
updated whenever a better solution is obtained. At termination 
all .arcs of Q have been examined and the incumbent gives the 

optimal solution. 
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However in case Q is constructed by priority order (c) # 
then the first time a incumbent greater than the upper bound 
of arc under consideration has been obtained gives the optimal 
solution. 

The procedure outlined in this section can be implemented 
by the following algorithm. 

2.51 Algorithm ADD'. 

* nrfc 'll** . ytit «u * «imw, v-*w 

INPUTS Network G (V, A), arc capacity c. . specific source s 

-t- J 

and sink t. Priority list Q ? upper bound on the arc 
potential of arcs of Q, Maxflow f . 

Stop Op. k 1 # f3c- 0 

Comment - ; f * denotes the current (maximum) additional flow 
obtained. 

Stop i s Select the k-th arc (i, j) from Q. IF f*> v^, GO TO 
Step 2, else 

IF f* < v. , then add ( i, j ) to G(V, A). 

Solve maxflow between s, t and let the flow be f^. 

g f k - £ 

IF g > f* then ft*. g and p-.- i, q^-D, G0 TO Step 2. 

Corrmien^s If priority set Q is according to the non-increasing 
order of upper bounds on arc potentials then at the 
first violation of condition > f , GO TO Step 3. 
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Step 2s k <■ k+1 

IF k > |Q| , GO TO Step 3, 

else GO TO Step 1. 

Step J5» (p, 0.) is the desired, optimal arc and the optimal 

increase in flow in the network G (V ? A) is f • STOP. 

2.52 Comp utational Perform ances 

The algorithm described above has been coded in FORTRAN- 10 
for DEC System 1090. A listing of the program is given in the 
Appendix. 

The program is tested on randomly generated problems. 

Test networks were generated as follows s The size par .Me tors 
i.e. number of nodes, number of arcs and length and breadth of 
the network, were fixed. Then a skeleton network consisting 
of arc disjoint path was constructed. The end points of the 
remaining arcs were generated randomly using uniform random 


number . 

The capacities of arcs J was generated using a uniform 

random number with range as 20 1 H(J) - T(J)i, 160 | H(J) - T(J)|, 

where H(J) and T(J) are the end points of the arc J. This per- 

^ „ .v,_ rnnaritv from uniform distribution 

mits the generation of the one capacity xrom ux 

with different ranges. 

,, <. 17 ^ rene-inm from 30 nodes to 200 nodes 

The problems of size ranging xium y 

i.l to .8 were generated. 


with density from 
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In Table 2.1 a comparison of problems of 30 nodes bo 
200 nodes with different arc set combinations is given. Alongwith 
the comparison of execution time for 3 priority rules the number 
of maxflow computations required for each priority rule is also 
given. For these problems the average execution time required 
by Pj is about 32 percent less than that required by 

As such we suggest that priority rule may be used. 


* 



CHAPTER III 


MULTIPLE ARC ADDITION IN A FLOW NETWORK 

3.1 INTRODUCTION 

In this chapter we will consider the problem of adding 
multiple arcs in the flow network with a specified source vertex 
s and sinc£ vertex t . 

We will like to select the arcs with minimum total capacity f 
whose addition will permit the flow to increase to the desired 
level. This problem is meaningful in the situations where the 
total cost of arc addition is a linear function of the total 
capacity of arcs. We shall further restrict ourselves to the 
selection of the arcs from the specified set E* which spans the 
min cardinality minimum cuts. Essentially it is the extension 
of the previous problems in which only one arc was permitted to be 
selected. 

We will consider a network G(V, A) where V „ s the node set 
and A is the arc set, with vertex Scathe source and vortex t as 
the sink. All arcs are undirected. Let c. . denote the capacity 

-L J 

of an arc (i, j), E' denotes the set of arcs which span the min 
cardinality minimum cuts. E' can be' obtained as described in 
Section 2.3* 
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A subset of E ! will be termed as feasible subset if the 
total sum of the capacity of the arcs in the subset is higher than 
the specified flow. 

A simple enumerative scheme for solving this problem can 
be constructed by solving max flow problem after adding each feasib] 
subset to the network. The optimal subset is a subset which increas< 
the flow to the desired value and for which total capacities of the 
arcs is minimum. The computational effort required for such a f 

scheme will be very high as there can be large number of feasible 
subsets of arcs in E' . 

We will describe an exact algorithm to solve this problem 
in Section 3.3. Exact algorithm is an implicit enumeration 
scheme and will require a procedure to obtain all subsets of a 
given set of elements. In Section 3.2 we describe one such 
procedure NEXKSB. This procedure is given in [14] . 

In Section 3.4 we will describe heuristic methods to solve this 
problem and will discuss computational performance in Section 3.5. 

3.2 PROCEDURE NEXKSB [l4]s 

Let us consider the combination of n things taken k at a 
time. There are (^) k subsets of an n-set. This procedure will 
generate them all in alphabetical order. In the lexicographical 
sequence, we obtain the successors of a given k-subset for n-set. 

For example let the n-set be Q = {l, 2, 3, 4,5}. The 3-subsets of 
Q in lexicographical sequence are given below? ■ 
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{1,2,3}, {1,2,41, {1,2,5}, {1,3,4}, {1,3,5}, {1,4,5}, 
{2,3,4}, {2,3,5}, {2,4,5} and {3,4,5}. 


In the lexicographical sequence, we- obtain the successor 
of a given k-subset {a^,..., a^} as follows. Search for the 
smallest h such that a k+ p_ k < n+l-h, then increase a k+ 2 _ k by 1 , 
and set — a j+i = k+ 2 -h,k). It is interesting to note 

that the index h can be found without searching. Indeed, at each 
transition from a set to its successors, h increases by 1 unless 
a k+l-h ^ n “k * ■'■ n wki 0 * 1 case h is reset to 1 on the next transi- 
tion. 


Algorithm NEXKSB [l4]s 

(1) m <s — 0; h« — kj GO TO 4 {first entry} 

(2) IF m > n-h, GO TO 3 ; h*-~0 {later entries} 

( 3 ) h<f — h +1 • m — a k+1 _ h 

(4) For a * 1, h: m+ l} 

If a 1 = n-k +1 GO TO 5 

(5) EXIT. 

3.3 EXACT METHOD FOR SELECTING ARCS FROM E’ s 

Next we shall describe an implicit enumeration scheme for 
selecting the minimum capacity arcs. 

Procedure IMSEARCH s 

oiW'HWMiairfflii rati. wi , 

Initializes Let the arcs of E 1 be ordered as 

^ a l , a 2* * * * » a m^ where 
c ( ) <. cG^) ^ ... < c( a m whcr -c 
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where cC?^) is the capacity of arc a^. 

For convenience we shall denote c(a^) hy c^ and the arc £3^ 
in this order by i. Further let g denote the 
increase in the flow desired. 

Step 12 Select, first k arcs such that when these arcs are 
added to the network, flow is Increased by atleast 
amount g . Let the capacity of these arcs be the 
current upper bound 

k 

UB = Z c . 

i=l 1 

For each subset P of the arcs(a^ } . . . , examine if 

C(P) + c k > g and 
C(P) + c k < UB. 

If yes solve a max flow problem and if the increase 
in flow p > g set UB = C(P). If all the subsets 
are enumerated go to Step 3 else repeat for next subset. 
Comments In the algorithm these subsets are generated such that 
first subsets of cardinality i = 1 are generated then 
subset of cardinality j + 1 and so on are generated. 

These subsets are generated using the procedure NEXKSB 
given in Section 3.2. 
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S tep 3 s k — k+l 

If k > m, go to Step 4 else select k-th edge from the 
list, add it to the list of selected edges and go to 
Step 2. 

Step 4? Stop. Output, list of selected subset, capacity, flow. 

In the next Section, we give a simple algorithm to 
execute this procedure. 

3,31 Algorithm IMSEARCHs 

INPUTS G( V, A) with c. . as capacity of arc (i,j), E' , g (required 

J-tJ 

flow increase) . 

Step Os Solve Maxflow (A, f) 

Comment: this procedure solves a maxflow problem between 
s and t in a network with arc in A and gives the maxflow 
as f . 

f «S f. 

S tep Is Arrange the arcs in B 1 in the nondecreasing order of the 

magnitude of their capacities. Let it be (a^, a 2 ».».,a n ). 
Let , E ~ (a^, , . . * , ^ 

R< — 0 , k M<—0, J <r-~l 9 C(R) 0. 

Step 2 s Select arc a k from E " . 

R^-~ R U {a k } 

C(R)-e— C(R) + c k 
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Step 3t If C(R) > g GOTO Step 4, Else k<~-k+l; GO TO Step 2. 

Comments C(R) is the total capacity of all arcs of set R. 

Step 4 s A U R 

Solve Maxflow (A, f) 
e „ — f - f * 

If e > g, NA~* k, U /c; C(R) 

GO TO Step 5A 

else ke^-k+1, GO TO Step 2. 

Step 5 A s Repeat for all i = 1, k - 1 
M-* — M U {a. } 

Step^5B, L <~--a k ^ A 0 a L 

If >_ g and < U 

Comments is the capacity of L-th arc in E 1 . 

Solve Maxflow (A, f), 
e< — f - f * 

IF e > g GO TO Step $6^ 

Else GO TO Step 6. 

Step 6 s Using Algorithm NEXKSB get next J- subset of M. Let it 
be X. If there is no subset of size J, GO TO Step 7* 

If c L + C(X) < U and > g 
n A< — A U X 

Solve Maxflow (A, f) 
e-fir--- f - f* 
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IF e > g, U«—c L + C(X), A« — A - X 
GO TO Step 6. 

S tep 7 « J — -J +1 

If J < NA GO TO Step 6, else GO TO Step 8. 

Step 8 s If }M| < n, k«s— k+1, «!«*=-• 1 

GO TO Step 5B, else GO TO Step 

Step 9& Outputs Current incumbent, U ,A, STOP. 

3.4 HEURISTIC METHODS $ 

In this section we will suggest the heuristics for solv- 
ing the problem given in Section 3.1. The first heuristic LBOND, 
essentially orders the arcs in order of their utilization, where 
utilization is the fraction of the capacity which is likely to be 
utilized when the arc is added to the network. The second 
heuristic UBOND is simpler in execution but may not give as good 
solutions as LBOND, 

H e ur istic J^BOND $ 

Initializes Let the network by G(V, A) . Required flow increase g t 
Maxflow in this network f* , E' the set of the arcs which can be 
added. 

Step 1? Compute lower bound p. . on arcs (i,j) £ S f using algorithm 
of section 2.49. 

Let F = { (i» j ) lp ± j = c^} 

order elements in F in nonincreasing order of c. . . 

» J".*] 
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R ^ — E - F 

c. . - p. . 

For (i,j) e R, compute d. . < — . 

J ±3 

Order arcs in R according to nondecreasing order of 
Let g be the remaining flow increase required and A 
be the arc set of the current network. 




Step 2 % If F =0 3 GO TO Step 3 else select an arc in F whose 
capacity is nearest to g' . Add this arc to current 
network arc set A. 

Solve Maxflow (A, f). Let e^-f - f *. If e > Cj GO TO 
Step 4, else remove the selected arcs from F, repeat 
Step 2. 

Step 3 % Select the first arc in R, add it to A, solve Maxflow 

(A, f). e 4r f - f* . If e > g, GO TO Step 4, 

else remove the selected arc from R, repeat Step 3 till 
R is empty. 

S tep 4 s STOP. OUTPUT, A, e 
3 . 41 Algor ithm LBOND i 

In this section we describe an algorithm for the heuristic 
in Section 3.4.- 
INPUTS G (V, A), c i5 , E r } } 

Solve Maxflow (A, f ) 

g <• — g 



38 


S tep I s Obtain lower bounds on mode potential on the network 
G(V, A) by method of Section 2.44. 

-1 

§Jj£IL-ir. s Repeat for all (i,j) e E' 

PlJ «— Hin {l i( 1 3 ) 

F - 0 
R< — 0 


Stepps For all (i, j ) such that (i, j) e E , 

If p ij ^ c id » 

then F<-F U {(i, j)}. 

Arrange all arcs in F in the nonincreasing order of 
their lower bounds on arc potentials. 

§te&it s For all (i, j) e E* such that ppj < c ± j. 




lid. ZFll 

c ii 


R^-R U {(i,d)} 


Arrange arcs in R in the non-decreasing order of 
capacity unutilization factor d. .. 

If F = 0, GO TO Step 7, else select the first arc (i,j) 
in F. If p, . >_ g' GO TO Step 6, else add (i, j) to A. 

X 

GO TO Step 8. 


S tep 5} 
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Step 6; 


— 1 + 1 

If ^ > | F | j GO TO Step 6B, else select the k’ -th arc,£. 


If p. . £ g' , repeat Step 6, else 

1 vJ 


i- i* 


St ep 6 b g Select t-th arc. Add it to A. GO TO Step 8, 

Step 7 » Select the first arc from R and add it to A. 

Step 8; Solve Maxflow (A, f) 
e — f - f* 

If e > g , GO TO Step 9, else remove (i,j) from E ' 


g'<-~-g - ' 

GO TO Step J. 




Step 9s Output (A). STOP. 

H euristic UBO NDs 

Initializes. Let E* be the set of arcs to be added, G (V, A) 

the current network, g as the required flow increase. 

A« — A E - E ' ) 

Step I s Compute upper bound on arcs of ET, and order them in 
order of nonincreasing upper bound. Let, 

E =-(&!» a 2 , a m }. 

Step 2 s Select first k arcs from E such that sum of their 

upper bound is greater than the required flow increase . 

Add these arcs to A. Let the new arc set be S’. Solve 

Maxflow (f , A) . Let e be the increase in the flow, then 

if e > ‘A , GO TO Step 3 , else remove these urcs from ET. 

k 

Let E = ff J U a, 

i«i k 

g V- • g “ e 
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GO TO Step 1. 

Step 5 s From A remove all those arcs of E 1 on which there is 
no flow . OUTPUT A. STOP. 
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3.5 COMPUTATIONAL PERFORMANCE i 

The exact algorithm and the heuristic LBOND described 
above has been coded in FORTRAN-10 for DEC System 1090. 

The program is tested on randomly generated problems. Test 
networks were generated as follows » The size parameters i.e. 
number of nodes, number of arcs and length and breadth of the 
network were fixed. Then a skeleton network consisting of all 
disjoint path was constructed. The end points of the remaining 
arcs were generated randomly using uniform random number. The 
capacities of arc J was generated using a uniform random number 
with range from 20 to 200. 

The problem of size ranging from 30 nodes to 150 nodes 
with density from 0.1 to 0.8 were generated. The cardinality of 
E was varied from very smaj.1 values ({E | = 5) to large values 
( |E ; | = 80) and for each problem at least two values of increase 
in flow were specified. 

In Table 3.1 a computational performance of these problems 
with different are set combination is given. The heuristic LBOND 
is fast compared to the exact method when the cardinality of E 
is large (> 20). However when cardinality of E is small the two 
methods perform indiscriminately. 

Also 94 percent of the solutions obtained by LBOND were 
within 25 percent of theoptimal value, and 62.5 percent of the 
solution obtained by LBOND were within 10 percent of the optimal 
value . 
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Table 3. is Computational performance of the exact and heuristic methods 



CHAPTER IV 


CONCLUSIONS AND RECOMMENDATIONS 

In this thesis we have considered the problem of expanding 
capacity of a flow network by adding arcs from a set of given 
arcs. We have considered two specific problems. In first 
problem, a single arc which increases the flow to maximum extend is 
required to be identified. ■ In the second problem, we select 

a { set of arcs from a given set of arcs such that the addition of 
these arcs can increase the flow capacity to a desired level and 
the sum of arc capacity is minimum among all such subsets. 

An implicit enumerative method has been developed for 
selecting single arc optimally to maximize the flow increase 
in an existing network. The computation of lower bounds on arc 
flows by the method of condensing the network gives useful infor- 
mation which can be utilized efficiently to select the order in 
which arcs are to be added in the enumerative scheme. Method can 
be further applied by obtaining better upper bound on arc flows 
which will help in reducing the number of maxflow problems to 
be solved to obtain the best arc. 

For the second problem, we have developed a implicit 
enumeration scheme and two heuristics. The implicit enumeration 
scheme for the multi arc addition problem gives efficient results 
in terms of computation time. Problem upto 150 nodes can 



be solved using this algorithm in less than 30 seconds. 

Since the problem considered in this thesis is a specific 
problem, a more general problem where arcs can be selected from 
any given set of arcs can be considered for further work. It 
may be noticed that the heuristic LBOND will work for any given 
set of arcs to be added, provided E / is replaced by the given 
set. For the exact method, however a new algorithm will have to 
be developed. 
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